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Abstract 

We use metastable NaCl-structure Tio.5Alo.5N alloys to probe effects of configurational disorder 
on adatom surface diffusion dynamics which control phase stability and nanostructural evolution 
during film growth. First-principles calculations were employed to obtain potential energy maps of 
Ti and Al adsorption on an ordered TiN(OOl) reference surface and a disordered Tio.5Alo.5N(001) 
solid-solution surface. The energetics of adatom migration on these surfaces are determined and 
compared in order to isolate effects of configurational disorder. The results show that alloy surface 
disorder dramatically reduces Ti adatom mobilities. Al adatoms, in sharp contrast, experience 
only small disorder- induced differences in migration dynamics. 
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Thin film growth is a complex physical phenomenon controlled by the interplay of ther- 
modynamics and kinetics. This complexity facilitates the synthesis of metastable phases, 
such as Tii_3.Ala.N alloys, which are not possible to obtain under equilibrium conditions and 
broaden the range of available physical properties in materials design. Fundamental under- 
standing of elementary growth processes, such as adatom diffusion, governing nanostructural 
and surface morphological evolution during thin film growth can only be developed by de- 
tailed studies of their dynamics at the atomic scale. Research has mostly been carried out 
using elemental metals, as reviewed in refs. [1,2]. Much less is known about the atomic-scale 
dynamics of compound surfaces, and particularly little about complex, configurationally dis- 
ordered, pscudobinary alloys which are presently replacing elemental and compound phases 
in several commercial applications. 

Kodambaka et al. [3] and Wall et al. [4, 5] used scanning tunneling microscopy to deter- 
mine surface diffusion activation energies Eg on both TiN(OOl) and TiN(lll). However, due 
to the vast difference between experimental and adatom hopping time scales, determining 
diffusion pathways requires theoretical approaches via first-principles methods that are capa- 
ble of providing clear atomistic representation on the ps time scale. Gall et al. [6] employed 
first-principles calculations to show that Eg for Ti adatom diffusion on TiN is much lower on 
the (001) than the (111) surface and used this diffusional anisotropy to explain the evolution 
of (111) preferred orientation during growth of essentially strain-free polycrystalline films. 
The correspondingly large differences in chemical potentials result in Ti adatoms having 
higher residence times on (111) than on (001) grains. 

Here, we use cubic Tii_^Al^N(001), a metastable NaCl-structure pscudobinary alloy, as 
a model system to probe the role of short-range disorder on cation diffusivities which control 
phase stability, surface morphology, and nanostructural evolution during growth. Tii-^Al^^N 
alloys with x ~ 0.5, synthesized by physical vapor deposition (PVD) far from thermodynamic 
equilibrium [7], are commercially important for high-temperature oxidation [8] and wear- 
resistant applications [9, 10]. Alloying TiN with AIN has also been shown to alter surface 
reaction pathways controlling film texture and nanostructurc [10-13]. Unfortunately, atomic- 
scale understanding of the growth of these important, and more complex, materials systems 
is presently rudimentary as best. Surface diffusion on a metal alloy, the CuSn system in 
ordered configurations and in the dilute limit [14], has only recently been considered using 
first-principles. However, it is well known that configurational disorder can have large effects 
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on the physical properties of sohd solutions [15] . 

We employ first-principles calculations using the projector augmented wave method [16] 
as implemented in the Vienna Ab-initio Simulation Package (VASP) [17], to determine 
the energetics of cation adsorption and diffusion on ordered TiN(OOl) and congurationally- 
disordered Tio.5Alo.5N (001) surfaces. Electronic exchange correlation effects are modeled 
using the generalized gradient approximation [18]. The plane wave energy cut-off is set to 
400 eV. We sample the Brillouin zone with a grid of 3 x 3 x 1 k-points. 

TiN(OOl), for reference, and Tio.5Alo.5N(001) surfaces are modeled using slabs with four 
layers of 3 x 3 in-plane conventional cells with 36 atoms per layer. Calculated equilibrium 
lattice parameters, ao, of bulk TiN, 4.255 A, and Tio.5Alo.5N, 4.179 A, arc employed. The 
vacuum layer above the surfaces corresponds to 5.5ao- The adatoms are spin polarized, 
which is found to be important for Ti adatoms with its partially filled 3d-shell, but not for 
Al. To investigate diffusion on a configurationally-disordered surface, the Tio.5Alo.5N(001) 
slab is modeled using the special quasirandom structure (SQS) method [19]. We impose a 
homogenous layer concentration profile and minimize the correlation functions on the first 
six nearest-neighbor shells for the slab as a whole. 

Convergence of diffusion barriers is tested with respect to the geometrical and numerical 
details of the calculations. Es results are within 0.04 eV of the converged value, partly due 
to error cancelation between the effects of treating Ti semicore states as core and the limited 
number of layers; both are of the order of 0.08 eV, but with opposite signs. Our primary 
focus is the observed differences in cation dynamics on the two surfaces. 

We begin by calculating the adsorption energy E^^f\x,y) for Ti and Al adatoms as 
a function of positions x and y on both ordered TiN(OOl) and disordered Tio.5Alo.5N(001) 
surfaces, 

E^n^, y) = Ef!;ll,{x, y) - - E^i^. (1) 

E'^ab+ad^^ the cucrgy of the slab with an adatom at {x, y), Egiab is the energy of the pure slab 
with no adatoms, and E^l^l is the energy of an isolated Al or Ti atom in vacuum. We use 
a fine grid of sampling points. Ax = Ay = 0.05ao. In each calculation, the adatom is fixed 
within the plane and relaxed out of plane. The upper two layers of the slab are fully relaxed, 
while the lower two layers are stationary. A periodic polynomial interpolation between the 
calculated points is used to obtain a smooth energy surface. 

Adsorption-energy profiles for Al and Ti atoms on TiN(OOl) and Tio.5Alo.5N(001) surfaces 
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FIG. 1. (Color online) Adsorption energy surface for: (a) an Al adatom on TiN(OOl), (b) an Al 
adatom on Tio.5Alo.5N(001), (c) a Ti adatom on TiN(OOl), and (d) a Ti adatom on Tio.5Alo.5N(001). 
Local minima are marked with red dots, while black dots indicate saddle point barrier positions. 
White lines, solid and dashed, marks preferred paths for diffusion on the disordered surfaces. 

are shown in Figs. l(a)-l(d). The most favorable sites for Al adatoms on both surfaces are 
directly above N atoms at bulk cation positions. For Al on TiN(OOl), Fig. 1(a), E^j^ is 
-2.54 eV. On Tio.5Alo.5N(001), Fig. 1(b), E^^ varies from -2.39 to -1.52 eV on the bulk 
cation sites depending on their local environment. Ti adatoms have two stable adsorption 
sites: fourfold hollows, surrounded by two N and two metal atoms, and the bulk site on-top 
N. For TiN(OOl), Fig. 1(c), E^l = -3.50 eV in the hollow site and -3.27 eV above N. On 
the alloy surface. Fig. 1(d), Ej^^ varies from -3.42 to -2.58 eV in the hollow sites and -3.23 
to -2.67 eV in on-top sites. Al-rich environments are much less favorable for both Al and 
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Ti adatoms as can be seen in the lower right regions of Figs. 1(b) and 1(d). The overall 
preferred sites for Ti on Tio.5Alo.5N(001) are fourfold hollow positions with one Ti and one 
Al nearest metal neighbors; not two Ti atoms as might have been expected. 

In order to quantify the impact of disorder on diffusion, transition state theory within a 
kinetic Monte Carlo approach is used to determine the mobilities of independent adatoms. 
The probability at each time step for a Ti or Al adatom at site i to jump to site j is calculated 
as 



where AEij — {Eij — Ej) is the difference between the adsorption energy in the local min- 
ima i and at the saddle point defining the barrier height Eij between sites i and j. The 
temperature T is 800 K, a representative value for PVD growth of transition-metal nitride 
thin films. For convenience, we choose the attempt frequency vq to be the same for Ti and 
Al on both TiN(OOl) and Tio.5Alo.5N(001) surfaces, but note that Al adatoms should have 
a slightly higher attempt frequency than Ti due to their lower mass. We determine nj(i), 
the probability density of finding adatoms on a given site i, at time corresponding to an 
ensemble average of a large number of individual cases. 

The most probable Al and Ti difi^usion paths are identified by imposing a constant prob- 
ability density of adatoms at the centers of circular grains with radii 8.5ao, and then prop- 
agating the probability density using Eq. 2. Adatoms crossing a grain boundary are not 
allowed to cross back. Thus, we obtain an adatom probability flow between sites i and j 
from the center of the grain outward. 



Steady-state results are plotted for Ti adatoms in Fig. 2 as grayscale intensity propor- 
tional to Fij. Panel 2(a) shows that the flow of Ti adatoms across the ordered TIN (001) 
surface is symmetric and utilizes all [110] paths. However, the flow of Ti atoms across 
Tio.5Alo.5N(001), panel 2(b), simulated using periodically repeated SQSs, is almost com- 
pletely absent in the energetically least favorable regions; most diffusion takes place along 
special paths. Such paths are indicated in Fig. 1(d) (Fig. 1(b) for Al adatoms) by white 
solid and dashed lines corresponding approximately to connections among the most favorable 
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local energy minima. 

Next, we determine the timescales of adatom diffusion on the two surfaces. Fig. 2(c) 
is a plot of the probability as a function of time that adatoms, individually placed at a 
randomly chosen site close to the center of a circular grain, have not yet reached the grain 
boundary. On the pure TiN(OOl) surface, Al and Ti adatoms show similar behavior as the 
somewhat higher barriers for Al diffusion are compensated by Ti adatoms having three times 
as many local minima positions for a constant grain size. The striking result, however, is 
that Ti adatoms diffuse much slower on Tio.5Alo.5N(001) than on TiN(OOl), while the rates 
for Al adatoms on the two surfaces are nearly equal. Since both Ti and Al adatoms diffuse 
predominantly along preferential paths on the disordered TiAlN(OOl) surface, the mobility 
differences are, in large part, explained by differences in energy profiles along these paths. 

Fig. 3 contains plots of E^^f^ relative to the most favorable adsorption site, along the 
preferred diffusion paths on Tio.5Alo.5N indicated in Figs. 1(b) and 1(d). The arrows in 
Fig. 1 define the starting position for the energy-path plots in Fig. 3. Corresponding 
Kds plots on TiN(OOl) are included for comparison. The calculated Al adatom diffusion 
activation energy on TiN(OOl) is Eg = 0.47 cV. Both the solid and dashed low-cncrgy paths 
for Al on Tio.5Alo.5N(001) exhibit the signature of configurational disorder with alternating 
deep and less-deep energy minima. However, the individual barrier heights are, in most 
cases, considerably lower on the disordered surface with the maximum barrier height just 
1.2 X larger than on TiN(OOl). E^ for Ti adatoms on TiN(OOl) is 0.40 eV and the smaller 
barrier for jumping out of the minima atop N is 0.17 eV. The individual barriers for Ti on 
Tio.5Alo.5N(001) are similar, but a series of less favorable energy minima, combined with 
asymmetric jump probabilities, creates additional migration obstacles approximately 2/3 
along the outlined paths. The maximum obstacles are 2. Ox and 1.6x higher than Eg on 
TiN(OOl) for the dashed and sohd diffusion paths, respectively, explaining the dramatic 
reduction of mobility in this case. The mass difference between Al and Ti atoms (which we 
ignored in these calculations) affects uq and will further increase the mobility difference. 

These results illustrate the complex effects that configurational disorder can induce on 
surface diffusion. They also help to understand aspects of the growth behavior of Tii-^rAl^rN 
thin films. Our observed increase in the residence time of Ti adatoms on Tio.5Alo.5N(001) 
vs. TiN(OOl) is consistent with the experimentally reported transition in texture for poly- 
crystalline TiN films, grown at relatively low temperatures with little or no ion irradiation, 
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from (111) [20] toward (001) upon alloying with AIN [10]. In addition, the higher mobility of 
Al, with respect to Ti, adatoms on Tii_^Al^N explains the results of Beckers et al. showing 
AIN enrichment in (111) and Al depletion in (001) oriented grains [11]. 

In conclusion, wc have compared the adsorption-energy landscape and the migration mo- 
bilities of Ti and Al adatoms on ordered TiN(OOl) and disordered Tii_^Ala;N(001) surfaces. 
The configurational disorder on the alloy surface results in the formation of deep trap sites 
for Ti adatoms which, together with an asymmetric adsorption energy map, dramatically 
decreases the Ti adatom mobility. In contrast, Al adatom mobilities are nearly the same on 
TiN(OOl) and disordered Tii_2.Ala.N(001) surfaces due to a much smaller disorder-induced 
spread in energy minima values and more symmetric diffusion probability distributions along 
the most favorable paths on the alloy surface. These results explain observed differences in 
preferred orientation and nanostuctural evolution during growth of polycrystalline TiN and 
Tii_3;Al3;N films. 
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FIG. 2. (Color online) Ti adatom diffusion paths from the center to the edge of (001) surfaces of 
(a) TiN, and (b) disordered Tio.5Alo.5N. (c) The probability as a function of time, that Al and Ti 
adatoms placed at random positions in the center of a circular grain of radius 8.5ao have not yet 
reached the grain boundary on TiN(OOl) and Tio.5Alo.5N(001). 
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FIG. 3. (Color online) Adsorption energies of Al (upper graph) and Ti (lower graph) adatoms 
along favorable diffusion paths on ordered TiN(OOl) and disordered Tio.5Alo.5N(001) surfaces. For 
the disordered alloy surface, energy profiles are plotted for both the solid and dashed paths across 
the SQS shown in Fig. 1. 
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